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ABSTRACT 

The outflow efficiency (if) from black hole (BH) accretion disc systems is known to de- 
pend upon both the BH spin (a) and the amount of large-scale magnetic flux threading the BH 
and disc. Semi-analytical flux-trapping models suggest retrograde BHs should trap much more 
large-scale magnetic flux near the BH leading to much higher r] than for prograde BHs. We 
self-consistently determine the amount of large-scale magnetic flux trapped by rapidly spin- 
ning (a = -0.9 and 0.9) BHs using global 3D time-dependent non-radiative general relativistic 
magnetohydrodynamic simulations of thick (h/r « 0.3-0.6) discs. We find that BH-trapped 
flux builds up until it is strong enough to disrupt the inner accretion disc. Contrary to prior 
flux-trapping models, which do not include the back-reaction of magnetic flux on the disc, our 
simulations show prograde BHs trap more magnetic flux, leading to about 3 times higher r] 
than retrograde BHs for \a\ - 0.9. Both spin orientations can produce highly efficient jets, 
rj ~ 100%, with increasing r] for increasing disc thickness. The similarity of r] for prograde 
and retrograde BHs makes it challenging to infer the sign of a based on jet energetics alone. 

Key words: black hole physics — (magnetohydrodynamics) MHD — accretion, accretion 
discs — galaxies: jets — gamma-rays: bursts — methods: numerical 



1 INTRODUCTION 

Various astrophysical systems — active galactic nuclei (AGN), bla- 
ck hole binaries (BHBs), and gamma-ray bursts (GRBs) — produce 
relativistic outflows. Mediated by large scale magnetic fields, out- 
flows can be powered by the central BHs via Blandford & Znajek 
(1977, BZ hereafter) mechanism or by the inner regions of accre- 
tion discs via Blandford & Payne (1982, BP hereafter) mechanism. 
Large-scale magnetic fields extract BH spin energy at a rate, 

PbZ= 4^°bh^h/(^h), (1) 

where k « 0.05 weakly depends on field geometry, ® B h = 
(1/2) jj e \B r \ 6A 6(p is an absolute magnetic flux through the BH, 
with the integral over all 6 and <p on the BH horizon, r H = r g {\ + 
Vl - a 2 ). Here r g = GM/c 2 , M, and a are BH gravitational radius, 
mass, and dimensionless spin, Q H = ac/2r H is BH angular fre- 
quency, /(0 H ) * 1 + L38(Q H r g /c) 2 - 9.2(Q H r g /c) 4 is a high-spin 
correction, dA 9(p = ^J II g d6d(p is an area element in 6 - <p plane, and 
g is the metric determinant (Tchekhovskoy et al. 2010, see also BZ; 
Lee et al. 2000; Komissarov 2001; Tanabe & Nagataki 2008). 

For fixed ® B h, Pbz is the same for prograde (a > 0) and retro- 
grade (a < 0) BHs. However, a massive accretion disc introduces a 
preferred spin/rotation direction. Does this lead to large differences 
in ®bh an d j et powers of prograde and retrograde BHs? Under- 

* E-mail: atchekho@princeton.edu (AT) 



standing this is crucial, as the link between jet power and BH spin 
has been confirmed in BHBs (Narayan & McClintock 2012, but see 
Fender et al. 2010) and jet power is increasingly often used for in- 
ferring BH spin (Daly 2011; Martinez- Sansigre & Rawlings 2011; 
Buliga et al. 2011; Lei & Zhang 2011; Bambi 2012). 

To quantify jet strength in an accretion system, we define BZ 
efficiency as jet power (eq. 1) in units of mass accretion rate, M, 

^ = 7^7 x 100% = ^0 2 H i^A /(0 H ) x 100%, (2) 
(M)c 2 47T \ c J 

where B h = ®Bn/((M)r 2 c) l/2 is BH dimensionless magnetic flux, 
and (...) is a time average. Similarly, we define magnetic flux en- 
closed by a toroidal ring, (r, 6), via <D(r, 6) = ff e , <0(p , B r dA &ip >, and 
its dimensionless version, 0(r,0) = 0(r,6>)/((M)^c) 1/2 , where the 
integral is over a polar cap, & < 6, of a sphere of radius r. 

Clearly, jet efficiency depends on the ability of accretion flows 
to drag large-scale magnetic fields to the centre. Lubow et al. (1994) 
studied magnetic field dragging by geometrically thin accretion 
discs, with angular thickness, h/r <§c 1. They found that field trans- 
port is governed by the value of the magnetic Prandtl number (Pr m , 
which is the ratio of turbulent viscosity to resistivity) and that in 
order for field dragging to be efficient, this number must be large, 
Pr m > (h/r)' 1 » 1. Magneto-rotational instability (MRI, Balbus 
& Hawley 1991) driven turbulent accretion flows have low values 
of magnetic Prandtl number, Pr m ~ 1 (Fromang & Stone 2009; 
Guan & Gammie 2009; Lesur & Longaretti 2009). For such natural 
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Retrograde BH, a = -0.9, 77 = 34% (model A-0.9f) 
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Figure 1. Initial conditions and time-evolution of our fiducial models, A-0.9f (left column) and A0.9f (right column). See Supporting Information for movies. 
[Panels a,f] A vertical slice through the initial conditions. Colour shows fluid frame rest-mass energy density (see colour bar) and thin black lines show levels 
of constancy of enclosed magnetic flux, (<D>, which represent field lines in the image plane. Magnetic field is axisymmetric, with = everywhere and 
B z > out to few x 100rg. Only the magnetic flux enclosed by the thick solid field line eventually falls into the BH. [Panels b,g] Show t- and (^-average of 
the magnetically-arrested state of the simulation. Solid lines show the same contours of <0> as in panels (a) and (f). Accretion accumulates so much flux in 
the centre that the inner disc is able to push only a fraction of the flux, enclosed by the thick line, into the BH. [Panels c-e,h-j] From top to bottom: rest-mass 
energy accretion rate, Mc 2 , dimensionless BH magnetic flux, 0bh, and energy outflow efficiency, rj. Both 0bh and 77 saturate ar t > 6000r g /c, beyond which 
the accretion flow is magnetically arrested. Dashed lines show time-averages: a prograde BH has a larger efficiency, 77 = 102%, than a retrograde BH, 77 = 34%. 



values, Pr m ~ 1, thin discs do not produce centrally-concentrated 
magnetic fields, which suggests that thin discs produce weak or no 
jets (Lubow et al. 1994; but see Spruit & Uzdensky 2005; Rothstein 
& Lovelace 2008). To include general-relativistic (GR) effects of 
the innermost stable circular orbit (ISCO) on field dragging in thin 
discs, Reynolds et al. (2006) proposed that gas and fields plunge 
into the BH inside the ISCO, so no magnetic flux passes through 
the "gap" between the BH horizon and the ISCO. Based on a GR- 
version of this flux-trapping "gap" model, formally applicable only 
to thin discs, Garofalo (2009); Garofalo et al. (2010) concluded that 
jets from retrograde BHs, a = -0.9, are > lOx more powerful than 
from prograde BHs, a = 0.9, both for thin and thick discs (see §2). 

Geometrically thick discs can efficiently drag large-scale 
fields inward even for Pr m ~ 1 (Cao 2011), and time-dependent 
simulations of McKinney & Gammie (2004); McKinney (2005); 
Hawley & Krolik (2006) show that thick accretion discs (h/r - 
0.2-0.3) can produce powerful jets. However, the simulated values 
of jet efficiency have a high degree of scatter, e.g., jet efficiency 
from retrograde BHs ranges in these works from 10% to 50% of 
the corresponding efficiency for prograde BHs at the same absolute 
value of spin. A major uncertainty in such studies is the depen- 
dence of jet efficiency on the value of large-scale vertical magnetic 
flux initially present in the disc, which is a free modeling parameter 
with no obvious natural value. Changes in the flux can significantly 
affect simulated jet efficiencies (McKinney & Gammie 2004; McK- 
inney 2005; Beckwith et al. 2008) and render them unreliable. 

Is it at all possible to obtain a well-defined value of jet ef- 
ficiency, free from the uncertainties in large-scale magnetic flux 



content of the initial simulation setup? Tchekhovskoy et al. (2011, 
TNM1 1 hereafter) showed that a promising approach is to start with 
a large vertical magnetic flux in the disc, more than the accreting 
gas can push into the BH. The excess flux remains outside, impedes 
the accretion, and leads to a magnetically-arrested disc (MAD, 
Bisnovatyi-Kogan & Ruzmaikin 1974; Narayan et al. 2003; Igu- 
menshchev et al. 2003; Igumenshchev 2008; TNM11). Estimates 
show that many astrophysical systems contain enough large-scale 
magnetic flux to naturally form MADs (Narayan et al. 2003; McK- 
inney et al. 2012). The inner disc properties of MADs are shown to 
be independent of the initial value of large-scale magnetic flux. So, 
we can reliably determine 77 for prograde and retrograde BHs and 
thereby test the flux-trapping "gap" model. In §2 we describe our 
numerical method and present our results and in §3 we conclude. 



2 NUMERICAL METHOD AND RESULTS 

We carried out time-dependent 3D general relativistic non-radiative 
MHD simulations using a high-order version of Godunov-type 
shock-capturing code HARM (Gammie et al. 2003; McKinney 
& Gammie 2004; Tchekhovskoy et al. 2007, 2009; McKinney & 
Blandford 2009) in modified spherical polar coordinates. We use 
logarithmically spaced radial grid, dr/r = constant, for r < r br 
(see Table 1 for r br values). For r > r br , the grid becomes progres- 
sively sparser, dr/r = 4(log r) 3/4 , with a smooth transition at r br . We 
choose grid inner radius, 7? in , such that there are at least 9 grid cells 
between the inner radial boundary and the BH horizon and place 
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Table 1. Simulation details 



Name a 


a 


7] b [%] 




A<p 


Resolution r 


in/fg 


r max l r g 


Rm/m 


Rout/ r g 


rbv/r g 


^run 


^avg i r g/ c ] 










Simulations with (h/r) 


rmax ~ 


0.2 and (h/r) 30 * 


n ^- 

U.J. 








A OQf 


KJ.y 


^zl + ^ 


1 nn 


2n 


988 v 1 98 v f\A 

ZOO A 1ZO A Kj'-t 


1 < 

1 J 


37.1 


7 


i n 5 

1U 


1 n^ 

1U 


(c\. 9nn^n^ 

(KJ, Z,KJKJOKJ J 


fsnnn- 9nn^»n^ 

^OUUU, Z.KJKJKJKJ) 


A HQ 
r\-\J.y 


n q 

— \J.y 


tO ± O 


1 nn 


7T 


988 v 1 98 v ^9 

ZOO A 1ZO A JZ 


1 c 

ID 


XI 1 

D 1 . 1 


n 8 

u.o 


in 5 

1U 


in 2 

1U 


(KJ, lOJ JD ) 


^snnn- 1 ^^^^^ 

^OUUU, lOJ JJ ) 


A n Q/ 


n q 

— \).y 


zK -a- 8 


1 nn 

1UU 




1zLzl v 198 v ^9 

Iti A 1ZO A jZ 


1 C 

1 J 


37.1 


8 


i n5 

1U 


1U 


(kj, 10 /Okj) 


r'snnn- 1 87^^ 

^OUUU, 10 / Ok)) 


r\-\J.7lQ 


— KJ.y 


ZL1 4- zL 


1 nn 


7T 


988 v f\A v ^9 

ZOO A 04 A JZ 


1 c 

ID 


XI 1 

D 1 . 1 


n 8 

u.o 


in 5 

1U 


in 2 

1U 


\KJ, 1 ID Ikj) 


^snnn- 1 i^icw 

^OUUU, 1 1 D Ikj) 


a n Q£i„ 

J\-\j.ynQ 


n q 


<4-U ± J 


1 nn 

1UU 


71 


988 v l^f\ v ^9 

ZOO A ZJO A jZ 


1 C 

1 J 


XI 1 

D 1 . 1 


n 8 
u.o 


i n5 

1U 


1U 


(\ 9^98- 1 
^IZjZO, lOJJJj 


(\ 9^98- 1 ^^^^^ 
^IZjZO, lOJJJj 


A n 

r\-\J.yH ( p 


— KJ.y 


jD ± J 


1 nn 


7T 


988 v 1 98 v f\A 

ZOO A 1ZO A Ot 


1 c 

ID 


XI 1 
D 1 . 1 


n 8 

u.o 


in 5 

1U 


in 2 

1U 


^snnn- aaia^ 

^OUUU, l^t / ) 


fsnnn- aaia 1 ^ 

i^OUUU, 1^1 / ) 


A H Qflin 
/\-U.>Tlip 


n q 

— KJ.y 


<4-U ± o 


1 nn 

1UU 




988 x 198 x ^9 


1 C 

1 D 


34.1 


0.8 


i n5 

1U 


m 2 

1U 


(\ Aim- 1 87Qn^ 

^l^-ZU/, 10 lyKJ) 


n Aim- 1 87Q^ 

^l"4-ZU / , 10 lyKJ) 


AH Of 




lUz ± 1 u 


1 nn 

1UU 


2n 


988 v 198 x 64 

Z.OO A 1 Z.O A U4 


1 c 
1 J 


34.1 


0.8 


in5 

1U 


m 2 

1U 


[KJ, lyo^J ) 


/'snnn- 1 Q8Q^^ 

(^oUUU, VyoyD ) 


AH Q 


n q 


Q£» -i-17 
yO ± 1 1 


1 nn 

1UU 




988 x 198 x ^9 

ZOO A 1ZO A JZ 


1 c 

1 D 


34.1 


0.8 


i n5 

1U 


m 2 

1U 


(f\. 999C<\ 
(kj, LLLoD ) 


fsnnn- 9998^"\ 

^OUUU, ZZZOJ ) 


AO Q7 


n q 

KJ.y 


113 + 13 


1 nn 




144 v 198 x ^9 

Itt A 1ZO A JZ 


1 s 


34.1 


0.57 


in 5 

1U 


10 2 

1U 


(c\- 9n^sn^ 

(KJ, Z.KJJDKJ) 


fsnnn- 9n^sn^ 

^OUUU, ZUjJU^ 


AO Q/i 


PI Q 


1 1 < _i_ 90 
1 1 J ± zu 


1 nn 


71 


576 x 128 x 32 


1 < 

1 J 


34.1 


0.8 


i n 5 

1U 


in 2 

1U 


(KJ, lOZOJ ) 


r'snnn- 1 f\if\^\ 

^ouuu, lOZOJ^ 


AO Qln 


n q 

KJ.y 


Q7 -i- 14 
y 1 zr. 1 


1 nn 


71 


288 x 64 x 32 


1 s 


34.1 


0.8 


in 5 

1U 


10 2 

1U 


(C\- 1 6Q1 

(KJ, iKjy ID ) 


/'snnn- 1 ^qi 

^OUUU, iKjylD) 


AO Q/^ 
rxKJ.yng 


PI Q 


1ZJ) ± zz 


1 nn 


71 


288 x 256 x 32 


1 < 

1 J 


34.1 


0.8 


i n 5 

1U 


m 2 

1U 


(\ A1C\1- 999Qn^ 
^14ZU / , LLLyKJ ) 


(\ A1C\1- 999Qn^ 
^IHZU / , LLLyKJ) 




0.9 


112+19 


100 


71 


288 x 256 x 64 


15 


34.1 


0.8 


10 5 


10 2 


(14207; 24410) 


(14207; 24410) 


AOMy 


0.9 


109 ±21 


100 


71 


288 x 128 x 64 


15 


34.1 


0.8 


10 5 


10 2 


(0; 15550) 


(8000; 15550) 


AQ.9h% 


0.9 


103 ± 14 


100 


71 


288 x 128 x 128 


15 


34.1 


0.8 


10 5 


10 2 


(8500; 14625) 


(8500; 14625) 


A0.9N25 


0.9 


116 ±32 


25 


71 


288 x 128 x 32 


15 


34.1 


0.8 


10 5 


10 2 


(0; 17350) 


(8000; 17350) 


A0.9N50 


0.9 


114 ±20 


50 


71 


288 x 128 x 32 


15 


34.1 


0.8 


10 5 


10 2 


(0; 14385) 


(8000; 14385) 


A0.9N200 


0.9 


118 ±32 


200 


71 


288 x 128 x 32 


15 


34.1 


0.8 


10 5 


10 2 


(0; 28620) 


(16000; 28620) 


A0.9R20 


0.9 


102 ± 19 


100 


71 


288 x 128 x 32 


20 


45.35 


0.8 


10 5 


10 2 


(0; 15295) 


(8000; 15295) 


A-0.94BfN30 
A0.94BfN30 


-0.9375 
0.9375 


88 ±75 
354 ± 43 


30 
30 


Simulations with (h/r) 
2n 136 x 64 x 128 
2n 136 x 64 x 128 


''max ~ 

10 
10 


0.6 and (/V r >30 ~ 
100 0.73 
100 0.73 


0.6: 
26000 
26000 


500 
500 


(0; 13000) 
(0; 13000) 


(8000; 13000) 
(8000; 13000) 



a Suffixes 6i hg" or "lg" indicate that ^-resolution was increased or, respectively, decreased by a factor of two as compared to the fiducial model. Models A0.9/#, 
A0.9hg, and A0.9/z#/^ are similar to model A0.9 but with 6- and (^-resolutions increased by factors fg and f<p, respectively, at t = I4201r g /c, where the pairs 
of values, (fg, f^), are given by (0.5, 1), (2, 1), and (2, 2), respectively. Model A-0.9flip is similar to model A0.9, but with BH spin value reversed to a = -0.9 
at t = I4201r g /c. Model A-0.9hg (A-0.9/^) is similar to model A-0.9 but with ^-resolution ((^-resolution) increased by a factor of two at t = 12328r^/c (at 
t = 8000rg/c, respectively). Model A0.9/zJ is similar to model A0.9/^ but with ^-resolution increased by a factor of two at t = 8500r^/c. 
b We quote 95.4% (two sigma) confidence intervals. c Given as NrXNgxN^, where Af r , Ng, and N<p are grid resolutions in r-, 6-, and (^-directions, respectively. 



the outer radial boundary at 7^ = I0 5 r g , which is larger than the 
light travel distance in a duration of the simulation (see Table 1 for 
R in values). This ensures that both radial boundaries are causally 
disconnected. We use standard boundary conditions (see TNM11): 
outflow in r-, reflecting in 6-, and periodic in ^-directions. 

We consider a retrograde fiducial model, A-0.9f, for spin 
a = -0.9, and a prograde fiducial model, A0.9f, for spin a = 0.9. 
Figure l(a),(f) shows our initial conditions: BH at the centre of an 
equilibrium hydrodynamic torus (Chakrabarti 1985; De Villiers & 
Hawley 2003), with angular velocity, Q oc r _1 75 . We place the torus 
inner edge at r in = 15r g and pressure maximum at r max ~ 35r g , 
so disc angular thickness at r max is (h/r) rmax « 0.2 (see Table 1). 

Here (h/r) r = [ff 6(p (0- 7i/2) 2 pdA e(p j ^pdAg^' 1 and the integral 
is over all ^ona sphere of radius r. We embed the torus with a 
purely poloidal (B^ = 0) loop of a weak magnetic field, with plasma 
P = Pgas/^mag > Anin = 100. The loop contains a large magnetic 
flux due to the large radial extent of the loop (see TNM1 1). 

We define rates of accretion of rest mass, M, and rest mass 
energy, F M = Mc 2 , by M = - ff e ^pu r dA e(p = F M /c 2 > 0, where 
p is fluid-frame mass density, u r is r-component of contravariant 
4-velocity, and the integration is over all 6, <p on the BH horizon, 
r = r H . Figure 1(c), (h) shows that after the start of the simula- 
tion M settles to a steady state at t > 5000r^/c. BH magnetic flux, 
0bh> increases until t ~ 6000r^/c, beyond which BH saturates (Fig- 
ure ld,i). Flux accumulates outside the BH, impedes the accretion, 
and leads to a magnetically-arrested accretion disc (MAD, Narayan 
et al. 2003). Some of the BH flux occasionally escapes from the BH 



via magnetic interchange (Stehle & Spruit 2001; TNM11), leading 
to oscillations of B h in time (Fig. ld,i). Figure 1(b), (g) shows the 
t- and ^-average of the flow over the MAD period (t avg in Table 1). 
In order to get to the BH, gas "diffuses" through the vertical mag- 
netic flux that reached size ~ 25 r g in a duration of our simulations. 
At large r > 30r g , all our discs have very similar values of angular 
thickness, (h/r) 30 « 0.3 (see Table 1). At smaller r, funnel mag- 
netic fields compress the disc vertically and lead to smaller h/r. 

We define total energy accretion rate (as measured at infinity), 
Fe = ff 0(p T r t dA 6ip , where T^ v is stress-energy tensor, the integral is 
over all 6, cp on BH horizon, and F E is positive if energy flows into 
the BH. We define energy outflow efficiency 77 as the energy return 
rate to infinity divided by the time-average mass accretion rate: 



(F M ) 

Time-dependence of 77 for our fiducial models, A-0.9f and A0.9f, 
is shown in Figure l(e),(j). The outflow efficiency saturates at 
t > 6000r g /c and clearly correlates with B h- This is not surpris- 
ing since we will see below that most of the outflowing energy is 
carried by BZ-driven jets, hence 77 « 77 BZ oc 0| r (eq. 2). A ret- 
rograde BH traps 30 units of BH (panel d), while a prograde 
BH traps 50 units (panel i). Correspondingly, the retrograde BH 
produces outflows (50/30) 2 « 3 times less efficiently (see eq. 2), 
with 77 = 34 + 3% (panel e), than the prograde BH, which has 
77 = 102 + 10% (panel j and Table 1). It is interesting to compare 
these results to those for thicker discs, (h/r) 30 « 0.6, and a some- 
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Figure 2. Equatorial profile of magnetic flux, (<f>(r, 6 = n/2)), in models 
different by the initial amount of large-scale magnetic flux (see legend). The 
ISCO radius is indicated by filled circles. Despite the magnetic flux content 
of the models varying by a factor of « 3, all prograde models (top cluster of 
curves) agree to < 10%. Hence, flux accumulation around the centre in our 
simulations is independent of the initial large-scale flux content of the disc. 



what larger magnitude of BH spin, \a\ = 0.9375. The retrograde 
model has rj = 88 + 75% and the prograde model has rj = 354+43% 
(see Table 1 and McKinney et al. 2012). Thicker discs reach higher 
77, and prograde BHs have higher 77 than retrograde BHs. 

We repeated our fiducial models in a reduced azimuthal 
wedge, A(p = n. We refer to them as A-0.9 and A0.9 and find 
77 = 43 + 8% and 77 = 96 + 17%, respectively, in agreement with 
our fiducial models. We carried out a number of resolution studies 
(see Table 1) and find they agree with our fiducial models within 
77 measurement uncertainty (see Table 1): the mean efficiency of 
our retrograde models is 77 = 38 + 4% and of prograde models is 
77 = 106+7%. Hence, 77 values are converged in our fiducial models, 
which at r = \0r g have Q z 95 and 20 cells per the fastest 
growing z- and <p- MRI wavelengths, respectively, and an equatorial 
cell aspect ratio, 6r : rd6 : rd<p « 2 : 1 : 7. Our highest (^-resolution 
model, A0.9/^, has Q z * * 100 and 6r : rS0 : r5<p * 2 : 1 : 2 
and is very well resolved according to MRI resolution criteria (e.g., 
Hawley et al. 201 1 ; Sorathia et al. 201 1). 

Why do outflow efficiencies of prograde and retrograde BHs 
differ? Can this be due to differences in initial conditions ampli- 
fied by accretion flow turbulence? To test this, we ran the prograde 
model A0.9 until t = 14207 'r g /c and instantaneously flipped the 
direction of BH spin. We refer to this model as A-0.9flip, which 
gives 77 = 40 + 8%, consistent with our previous retrograde results. 
But what about large-scale magnetic flux? Since its value is con- 
served, the differences in the flux cannot be erased by turbulence; 
can they affect the outcome? To test this, we carried out a series of 
simulations different only by the magnitude of the magnetic field, 
i.e., with p m \ n = {25, 50, 100, 200}. We refer to these simulations as 
models A0.9N25, A0.9N50, A0.9f (our fiducial prograde model), 
A0.9N200, respectively. The initial fluxes in these models differ 
by a factor of « 3. Does this lead to a similar difference in the 
magnetic flux that reaches the BH? Figure 2 shows that all models 
have equatorial magnetic flux profiles, (0(r, 6 = tt/2)), that agree to 
< 10%, and 77 that agree to < 20% (see Table 1). This demonstrates 
that 0bh and 77 are independent of the initial magnetic flux content 
of the flow. We also verified that our results are insensitive to the 
initial position of the torus (see model A0.9R20 in Table 1). 

So far we considered the total energy output of the BH. What 
fraction of this energy goes into jets and into winds? The t- and (p- 
averages of our fiducial models, A-0.9f and A0.9f, shown in Fig- 
ure 3, have similar hourglass shapes. The equatorial part of the ac- 
cretion flow directly reaches the BH, while the inflow at higher lati- 



tudes turns around and forms weakly magnetized disc winds. These 
winds confine the highly magnetized relativistic polar jets that di- 
rectly connect to the BH. Figure 3(a) shows that in the retrograde 
model, A-0.9f, the efficiency of the wind is 77 wind (-0.9) « 7%. The 
remainder comprises jet efficiency, 77j et (-0.9) « 34 - 7 = 27%. 
Similarly, Figure 3(b) shows that wind efficiency for the prograde 
model, A0.9f, is 77 wind (0.9) « 16%, so jet efficiency is //jet (0.9) « 
102 - 16 = 86%. In both cases, the jet, which is predominantly 
powered by the BZ mechanism, carries most (^ 80%) of the total 
power output, and the power of the wind, launched by a BP-like 
mechanism, is subdominant. Note that all energy flow streamlines 
start on the BH and none start in the body of the disc. This is a man- 
ifestation of energy conservation: energy can either come from the 
BH via the release of gravitational binding energy or Penrose/BZ 
effect, or it can be advected inward from large distances by accret- 
ing material (if the accretion flow is unbound at large radii). 

Our results differ from semi-analytical flux-trapping "gap" 
models (Garofalo 2009), which predict for \a\ = 0.9 that retrograde 
BHs produce > lOx more powerful jets than prograde BHs. The 
"gap" models make use of a sharp transition at the ISCO to a plunge 
flow that sweeps into the BH all flux from an area (the "gap") be- 
tween the BH horizon and the ISCO. Since retrograde BHs have 
the largest ISCO area, in the "gap" model they receive the largest 
flux and produce the most powerful jets. However, poloidal flux 
distributions in our prograde and retrograde simulations are strik- 
ingly similar and show no "gap" (see Figures 2, 3): clearly, the 
"gap" model does not apply to thick radiatively-inefficient discs 
(unlike used by Garofalo 2009; Garofalo et al. 2010). This is not 
surprising because the sharp transition to the plunge at the ISCO, 
which the "gap" model relies on, is pronounced only for very thin 
(h/r < 0.05) radiatively-efficient discs (Penna et al. 2010). In ad- 
dition, a key piece of physics is missing from flux-trapping "gap" 
models: "gap" models assume that magnetic field is weak and has 
no back-reaction on the plunging inflow. However, BH magnetic 
flux builds up to a natural saturation point at which the field does 
back-react: the magnetic flux not only modifies the rotation rate of 
plasma in the plunging region, it also escapes from the BH via mag- 
netic interchange (Stehle & Spruit 2001; TNM11; McKinney et al. 
2012). In our simulations, this non-axisymmetric mechanism fills- 
in the ISCO region with magnetic flux. No "gap" forms because the 
inflow of magnetic flux happens just as fast as the outflow of mag- 
netic flux, unlike in the flux-trapping "gap" models. "Gap" models 
have not yet accounted for these effects, which dominate in MADs 
and, hence, in accretion systems with the most efficient jets. 



3 CONCLUSIONS 

We set out to determine whether prograde or retrograde BHs 
produce the most efficient outflows. Accretion discs with large- 
scale vertical magnetic fields naturally evolve into a magnetically- 
arrested disc (MAD) state in which the central BH is saturated with 
magnetic flux and the outflow efficiency (77) is maximum. Our simu- 
lations show that in this state, 77 only depends on the BH spin, a, and 
the angular thickness of the accretion disc, h/r, and is independent 
of the initial magnetic flux content of the disc (see §2). 

We find that for medium-thickness discs, h/r « 0.3, and an 
absolute value of spin, \a\ = 0.9, a retrograde BH produces jets and 
winds with an efficiency, 77 = 38 + 4%, which is a few times smaller 
than a prograde BH, 77 = 106 + 7%. In both cases, most of the 
energy (^ 80%) emerges in the form of BZ-powered, relativistic, 
highly magnetized collimated polar jets (see §2). We expect the 
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Retrograde BH, a = -0.9, t?=34% (model A-0.9f) Prograde BH, a=0.9, r? = 102% (model A0.9f) 




Figure 3. Time and azimuthal average flow pattern in our fiducial models, A-0.9f (panel a) and A0.9f (panel b). BH is shown with black filled circle. Black 
and white dashed lines show poloidal magnetic field lines (constancy levels of dimensionless magnetic flux, = {10, 20, 30, ...}). Gray shading shows velocity 
streamlines. Green lines show the position of stagnation surface, at which u r = 0. Colour shows energy outflow efficiency enclosed between the point in 
question and an equatorial energy flow streamline (see colour bar), so the edges of colour bands represent energy flow streamlines. The flow pattern is a 
standard hourglass shape: equatorial inflow at low latitudes, inflow-outflow at intermediate latitudes, and twin polar jets (shown in red) at high latitudes. 



above values of rj to serve as upper bounds on jet efficiency for 
lower BH spins, \a\ < 0.9, and thinner discs, h/r < 0.3. We also 
find that a two-fold increase in disc thickness (to h/r « 0.6) leads 
to about a three-fold increase in rj (see §2). These results can be 
used to place limits on the spin of central BHs from the observed 
values of jet efficiency. However, based on jet energetics alone it 
will be challenging to tell apart prograde and retrograde BHs due 
to their similar efficiencies. Future studies should consider radiative 
effects, which are not included in the current work, and investigate 
spin-dependence of 77 for thinner discs, with h/r < 0.1. 
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SUPPORTING INFORMATION 

Additional Supporting Information may be found in the online 
version of this article: 

Movie files. Movies of models A-0.9f and A0.9f (click for movies). 
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